Accounting for the Unresolved X-ray 
Background with Sterile Neutrino Dark Matter 



D.T. Cumberbatch and Joseph Silk 

Department of Astrophysics, University of Oxford, Keble Road, Oxford, 0X1 3RH 

Abstract. We consider a scenario where keV sterile neutrinos constitute all of the currently inferred 
dark matter abundance, whose radiative decays could potentially account for the flux contributions 
to the X-ray background (XRB) by unresolved sources. Here we apply integrated flux methods to 
results from the observations of the North/South Chandra deep fields (CDF-N/S) in order to deduce 
constraints on the sterile neutrino mass-mixing parameters. 
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INTRODUCTION 

The proportion of the cosmological energy density manifested as dark matter is now 
well determined by observations to within 10% [yj], yet its nature remains a subject 
of extensive debate. Recent astrophysical observations, such as indications of central 
cores in low-mass galaxies (see e.g. [0]), and the deficiency of satellites observed in 
Milky Way-sized galaxies (see e.g. [3D), indicate possible shortcomings of the [A]CDM 
paradigm and have boosted interest in a warm dark matter (WDM) scenario which may 
alleviate these "small scale problems" (see e.g. [4]). One of the most popular WDM 
candidates are right-handed ( "sterile" ) neutrinos v s , since they naturally arise in many 
extensions of the Standard Model (see e.g. j^). Other motivations for sterile neutrinos 
to possess mass-mixing parameters necessary for them to be produced via oscillations 
(with the familiar "active" neutrino species) include the enhanced production of neutral 
hydrogen before reionization [6] and a possible explanation for the baryon asymmetry 
of the Universe [5Q. 

Sterile neutrinos can be produced through non-resonant oscillations, the simplest 
being the popular Dodelson-Widrow (DW) mechanism [7], or resonant oscillations 
for cosmologies which include a finite lepton number asymmetry. In both scenarios, 
sterile neutrinos are predominantly produced during the QCD epoch. Unfortunately, 
uncertainties relating to hadronic interactions during this time result in uncertainties in 
the relationship between the mass-mixing parameters and the relic abundance of sterile 
neutrino dark matter 

Direct constraints on the mass-mixing parameters can be obtained by exploiting X-ray 
observations since keV-mass sterile neutrinos posses a radiative decay channel, yielding 
photons potentially detectable in astrophysical X-ray sources, including the diffuse 
cosmic X-ray background [9], the Milky Way liol fllll as well as being inferred from the 
X-ray spectra of nearby galaxies [12] or clusters (see e.g. 10]). Indirect constraints can 
be achieved from Ly-oc forest measurements II 1411 . since they are sensitive tracers of the 



primordial density fluctuations on the smallest scales where WDM typically suppresses 
clustering. 

Assuming that sterile neutrinos constitute all of the dark matter and are produced 
solely through the DW mechanism, recent X-ray analyses (see e.g. [@, Qjl Qli E2L 
T3IP deduce upper limits on their mass, m s ~ [3 — 8] keV. Conversely, the latest Ly-a 
analyses [15] of the high redshift power spectra from the SDSS iflill furnish constraints 
m s < (10—13) keV. Hence, these results seem to exclude scenarios where all dark matter 
consists of sterile neutrinos produced through the DW mechanism. (However, such 
evidence does not preclude the possibility of a sub-dominant sterile neutrino dark matter 
component [1J/ZD •) However, hidden systematic effects within SDSS measurements of 
the Ly-a flux power spectrum could relax the severity of the corresponding constraints 
inferred from the 1-D matter power spectrum from the Ly-a forest 111 in . 

In these proceedings we apply an integrated flux method to the CDF-N/S spectra of 
the XRB in order to constrain the radiative decay rate of DW sterile neutrino dark matter. 



FLUX FROM STERILE NEUTRINO DECAY 

Sterile neutrinos possess a radiative decay channel into an active neutrino and a photon 
with energy Ey = m s /2. For Majorana neutrinos, the radiative decay rate IJilQji] can be 
expressed as 

r y ~1.38xlO~ 22 sin 2 20,f T ^ 7 ) s -1 , (1) 



IkeV/ 

where 6 S is the mixing angle in vacuum between the sterile and the active neutrino 
species. (For Dirac neutrinos the above rate is halved 

The predicted signal from such decays originates from both extra-galactic (EG) neu- 
trinos, as well as those residing in the Milky Way (MW) halo. The EG contribution can 
be evaluated assuming a uniform distribution of neutrinos in the visible Universe up to 
very small redshifts. The differential energy flux (energy flux per unit energy, integrated 
over the field of view (FOV), ^fov, of the detector ) can be expressed as []23j] 

g EG ^ r 7 Qdm Pc QFOV ^) 

9e ~ A%m s H{m s /2E 7 -\) ' 

where p c is the present critical density, Qdm — 0.21 is the dark matter density parameter, 
and H(z) is the Hubble function, where in this study we adopt a flat A- matter dominated 
universe with H(0) = 73 km s _1 Mpc -1 , as well as matter and dark energy density 
parameters equal to Q m ~ 0.24 and Q./^ ~ 0.76 respectively yj]. 

The sterile neutrinos in the Galactic halo give rise to the differential energy flux 



pMw = _y_ j j ( Q)dxdQ = _L_ f 5dm (a) d a = ^a FOV 5dm, 

Anm s Jfov Acs. 4xm s Aov 4xm s 

(3) 

which depends on the integral over the dark matter density pd m along the line of sight 
(l.o.s.), 5dm.- Owing to the small variation of 5d m over the small FOV of Chandra (~ 5' 
circular area) in the directions of interest, we can approximately replace Sd m (i-e. the 



average value of Sdm over the relevant FOV), with the value of S am at the centre of each 
FOV. 

We calculate the flux contributions from decays within the Galactic halo using the 
Navarro-Frenk-White (NFW) density profile p d m(r) = p, (r/r s y l (1 + (r/r s )y 2 
where p s and r s are the scale density and radius respectively. We adopt the recent 
evaluation of the virial mass of the MW halo, 0.6 x 10 12 M Q < M vir . < 2.0 x 10 12 AT©, 
from [|22|l . which is also consistent with i? v i r . = 255 kpc, R w i T ./r s = 18 and R & = 8 kpc. 

Adopting the above, we obtained values for the dark matter surface density of 
(Sdm)cDF-N = [0.0115,0.0384] gcm~ 2 in the direction of CDF-N (with Galactic coor- 
dinates (Z,fc)cDF-N = (125.89°, 54.83°)), and (S dm )cDF-s = [0.0111,0.0369] gcm~ 2 in 
the direction of CDF-S (with Galactic coordinates (/,Z?) C df-s = (223.57°, -54.44°)), 
where the lower (higher) extremes of each the ranges stated correspond to the lowest 
(highest) values of the MW virial mass mentioned above. 



RESULTS 
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FIGURE 1. Left - XRB spectra in terms of the total specific intensity I v = dq>E / dQ. (thick-black crosses), 
and contributions from unresolved sources (thin-black crosses (CDF-N), thick-grey crosses (CDF-S) ). 
The errors displayed are taken to be at the 1(7 level. Right - 1 <y upper limits on sin 2 (20 4 ) using values 
of the Galactic halo virial mass of M v ; r . = 0.6 x 1O 12 M (dashed lines) and M v ; r . = 0.6 x 10 12 M Q (solid 
lines), for CDF-N (black lines) and CDF-S (grey lines). The dot- dashed line indicates the parameter space 
permitted for DW sterile neutrinos constituting all of the currently inferred dark matter relic abundance. 

Here, we utilise results from the analysis of the CDF observations by Worsley et 
al. 112 111 . Fig.[T](/e/iO displays the total XRB spectra aswell as the contributions attributed 
to unresolved sources. The errors displayed are taken to be at the la level. We then 
derive upper limits on the radiative decay rate of sterile neutrinos by applying the 
following spectral analysis. 

For each selected value of m s , we add the corresponding EG and MW contributions 
(after accounting for signal broadening, owing to Chandra's finite energy resolution, 
(which we conservatively estimate as AE/E > 0.03 [24]), and limit the decay rate (0Q) 
by invoking the following criterion 



f f v DF - N/S dE> I (l™+I™ w )dE 

JAE dp JAE dp 



(4) 



where the left-hand integral is evaluated over the upper limit of each of the data points 
displayed, each with energy range AE^ p . In fig. \T\ (right), we display the resulting la 
upper limits on sin 2 (20s) as a function of m s , obtained using values of the Galactic 
halo virial mass M v i r . = 0.6 x 10 12 M Q (dashed lines) and M v i r . = 0.6 x 1O 12 M (solid 
lines). In addition, we plot a contour indicating the parameter space permitted for sterile 
neutrinos constituting all of the currently inferred dark matter relic abundance when 
generated solely via the DW mechanism (dot-dashed line). 

Hence, using fig. Q] (right), the la upper mass limits for DW sterile neutrino dark 
matter are 

f [5.95,8.30] keV for CDF- N observations 
Ms < \ [6.65,9.49] keV for CDF-S observations 

where once again, the lowest (highest) extremes of each mass range correspond to the 
lowest (highest) values mentioned above for the virial mass of the Galactic halo. 

In conclusion, the above constraints are slightly less stringent than those obtained 
from other related works which utilise a line signal non-detection spectral analysis of 
the XRB data (see e.g. llllll ). However, we consider that such techniques yield overly 
stringent results when the energy resolution of the data is not sufficient enough to 
resolve line signals from local contributions. If so, then such results should be viewed 
as absolute lower/upper bounds on the sterile neutrino mass/mixing-angle, whereas the 
results obtained using the above integrated flux method should be considered to be the 
converse bounds on these parameters. 
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